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S : ABSTRACT 

We discuss the weak coupling expansion of massless QCD with the Dirac operator which 
is derived by Neuberger based on the overlap formalism and satisfies the Ginsparg- Wilson 
relation. The axial U{1) anomaly associated to the chiral transformation proposed by 
Liischer is calculated as an application and is shown to have the correct form of the 
topological charge density for perturbative backgrounds. The coefficient of the anomaly is 
evaluated as a winding number related to a certain five- dimensional fermion propagator. 
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Recently Dirac operators which may describe exactly massless fermion on a lattice 
are proposed by Neuberger |I| based on the overlap formalism ^ and by Hasenfratz, 
Laliena and Niedermayer |Q based on the renormalization group method These 
Dirac operators are very different from each other but both satisfy the Ginsparg- Wilson 
relation [0, which ensures that the fermion propagator anti-commutes with 75 at non- 
zero distances and thus respects chiral symmetry in that sense. Subsequently Liischer 
observed that the Ginsparg- Wilson relation implies an exact symmetry of the fermion 
action and the anomalous behavior of the fermion partition function under its flavor- 
singlet transformation is expressed in terms of the index of the Dirac operator arised as 
the Jacobian factor of the path integral measure, providing a clear understanding of the 
exact index theorem on a lattice in Ref. 

Though it was pointed out in Ref. that Dirac operators satisfying the Ginsparg- 
Wilson relation will necessarily exhibit non-local behavior in the presence of the dy- 
namical gauge fields, a recent perturbative analysis of chiral gauge theory in the overlap 
formalism |0 suggests some of them will be well controllable as a field theory on a lattice. 

In this paper we consider the weak coupling expansion of massless QCD with the 
Dirac operator which is derived by Neuberger and satisfies the Ginsparg- Wilson relation. 
Then we apply our results to analyze the axial U{1) anomaly following the formulation 
of Ref. and confirm that the Jacobian factor is given in the form of the topological 
charge density for slowly varing perturbative gauge fields, supplementing the original 
calculation of the same quantity of Ref. |^ . 

The covariant anomaly in the vacuum overlap formalism has been discussed by Neu- 
berger and Narayanan Randjbar-Daemi and Strathdee |jlO| and Neuberger 



Note also that the axial anomaly has been calculated in the context of the domain-wall 



QCD by Shamir p| 



We begin with a brief review of the symmetry argument of Ref. The partition 
function we consider is of the form 

"##6"'^'^'^^^ (1) 
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with a lattice Dirac operator D satisfying the Ginsparg- Wilson relation 



D75 + ToD = aD-f^D. (2) 

In Ref. [|] Liischer pointed out that with the aid of the relation the fermion action 
a^TjipDip is invariant under the infinitesimal transformation ip ^ ip + ieTSip, ip 
ip + iedipT where T is a generator of the rotation in the flavor space and 

5^ = 75(1 - ^aD)7/;, 5^ = i,{l-hxD)^^. (3) 

The path integral measure yields the Jacobian factor —ieati{T'-)r^D}, which does not 
vanish only for the flavor singlet chiral rotation, accounting for the index theorem of 
Ref. i. 

A Dirac operator satisfying the relation was proposed by Neuberger in Ref. [|l| 
and now we review its brief derivation starting from the domain wall fermion of the 



Shamir type ||T3| based on Ref. [|T^], trying to clarify the physics backgrounds of the 
Dirac operator. 

The domain wall fermion of this type is a Wilson fermion in five dimensions with the 
finite size of the fifth space and is described by the action 

S = ^a^X!^("^''5)^5(Ai)ms,n<^(n,t), (4) 
s,t rn,n 

^ 1 1 

DMms,nt = 7mCm("^' '^)^s,t + Br,{m, n)Ss,t ^L<5s+i,t PrSsmi 

+HPLSs,NSt,l + fJ'PRSs,lSt^N, (5) 



1 



Cf,{m,n) = — 6^+f,^nUf,{m) - 6m,n+f,Ul{n) , B5{m,n) = — + B{m,n), (6) 
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as 



B{m,n) = — + [^^m,n - Sm+i,,nU^{m) - 5m,n+f,Ul{n) , (7) 

a , , 

where m, n and a denote the four dimensional space indices and their lattice spacing 
while s,t and as denote the indices and the lattice spacing of the fifth dimension (1 < 
s,t < Nr,). In the action Eq. the link variables f/^(m) couple only between the fields 
iplm, s) and ilj{m ± /i, s) so that there is no gauge interaction in the fifth space. The 
terms proportional to 1/as come from the kinetic and Wilson terms in the fifth direction 
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and the Wilson parameter in that direction is set to be unity. We have omitted the 
factor as in front of the sum J2st so that ip has the correct dimension 3/2 for fermions 
in the four dimensions. The role of the parameter /i is explained later. 

The action Eq. (^) is also regarded as the action of the A^s flavor Wilson fermions in 
four dimensions with a specific mass matrix. With a suitable choice of Mq the action 
describes one Dirac fermion with the mass /i + (9(e~^^/a) and A^s ~ 1 Dirac fermions with 
the mass of the order of the inverse lattice spacing |T5|. Then the following partition 
function 



= det L)s(/i = 0)/ det L)s(yU = — ) (8) 

as 

may describe one massless Dirac fermion regulated a la Pauli-Villars by one Dirac 
fermion with the mass l/as if the A^s ~ 1 heavy fermions cancel out. Such cancellation 
allows to take the limit oo. Then taking the limit as — > 0, one massless Dirac 

fermion remains. 

Now we compute the determinant det-Ds(Ai) and obtain the final expression of Eq. 
(§) as a single determinant, which allows the path integral expression like eq. In 
the chiral basis of the 7-matrices defined by 

7m = ) ' 75 = J ^ , f^M = (-^^1' -^^2, -«CT3, 1), (9) 

the Dirac field is written in the chiral components as 

V'=(^^), ^ = (v^h,V'l). (10) 

We also introduce the notation 

= (11) 
althoug it does not mean hermitian conjugate. Then the action can be written as 



B, \ ( \ 
B, -C r^'+[ -l/as r'''^' 



( -l/as \. ./'oMaaj./'O'^^aaI/' Mn,t) 



(12) 
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Thus 'j4^D5{fi) takes the form 
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(13) 



where = J2ii'^l.Cfj,{'m,n) and — C = ^^<y^C^{m^n). Assuming the four dimensional 
space consists of sites and the color and flavor groups are SU{Nc) and SU{Nf ), each 
block C, -B5 etc. above is g x g matrix with q = 2NcNfL^. The determinant of this 



matrix is evaluated following the technique developed in Ref. [0, 0. Moving the first 
q columns to the last, 74-D5 takes the from 
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(14) 



Decomposing the matrix (0) into the product of the two matrices as 
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f3i \ 
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vx = -ai ^Pi 
-PiVi_i - ttiVi = 



(15) 
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its determinant is given by 

det ttj det(l — f ATg) (16) 
and the relations at the end of eq. (|15]) leads to 

vns = (-"^s/^^s) ■ ■ ■ (-ar^A)- (17) 

In the case of our 74-D5, 

and for the later convenience we introduce the transfer matrix T and the Hamiltonian 

as 



a 13-^,1 74, J - e ' ^ " |^ B^'C^ a,C^B,'C + a,B, 
The final expression is 

det^,D,{fi) = (-l)«(^-i)(a5)-''^^(deti?5)^^ 



(19) 



(20) 



where (— 1)'J(^5 arises when the first q columns are moved into the last column and 
eq. (P) is given by 



det ^{1-75 ta^nhi^N^a^H^)}. (21) 



Taking the limit A^^s ^ 00, 



tanh(^iV5a5i/5) ^ e{a,H,) = . (22) 



In the limit 05 — > 0, 



where X is the Wilson-Dirac operator on a lattice. Therefore in this hmit, eq. is (up 
to constant) 

detA /^ = -(l-75^) (24) 

which allows the path integral expression over the fermion fields as eq. (|l|). The 
Ginsparg- Wilson relation is reduced to {H/\fW^Y = 1 for Dirac operators of the form 
eq. (|2^) , which is in fact satisfied by definition fl^ . 

Now we discuss the weak coupling expansion of the Dirac operator: 

aD = 1 + X , ^ :, (25) 

Xnm = i^Cf,{n,m) + B{n,m) - -M^dnrn- (26) 
The Wilson-Dirac operator X 

may be expanded as 

X(g,p) = Xo(p) (27r)V(g -p) +Xi(g,p) +X2(g,p) + 0{g^), (28) 



where 



^oip) = -7m sin ap^ + - ^ (1 - cos ap^) Mq, (29) 



a CL~~ a 



Xi{q,p) = J-0^,{27rr6{q-p-k)gA,{k)V,,(p+'iy (30) 

^A,(A:0A,(fc2)l^2,(p+^). (31) 
The vertex functions are given explicitly as 

Vi^i + I) = ^^/^ a (^Pt^ + + r sin a |^p^ + (32) 

-^0 b + 



dp, " V 2; ' 

V2^,(p+^Y'j = -i7^a sin a + yj + ar cos a + yj . (33) 



Let us assume that 1/VXtX can be expanded in the following form, 

1 . . / 1 ~ 



Then it should satisfy 



{p, q) + Yi{p, q) + 1^2 (P, q) + 



(XtX)(p,g) 
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is,t) = Sip-t). 



(34) 



(35) 



We have made use of the abbreviations as / = Jf^ and 6{q — p) = (27r)^(5^(g — p). In 



(s,t) 



the first non-trivial order, the above equation reads 

i (v^„(p. t) + ^ /(.f t.Y)„(p, ,)F,(,, .) (-^ 

+ /(A'tA')^(p,,)(^)„fe.i)=0. (36) 

Since {y X'^X^ ^ {p, q) is completely diagonal with respect to the momentum, spinor and 
color indices, 

(yxtx)jp,g)=a;(p)5(p-g), (37) 

where 

auij)) = ^ sin^ap^ + ^^^(-'^ ~ cosap^) — Moj > 0, (38) 

Eq. (^) is reduced to the relation between the elements Yi{p,t) and (X"l'X)i(p, t), 
leading to the solution |[T7[| , 

i'i(p.«) = 



u;(p)a;(i:) '^ci;(p) + ci;(t) ■ 
Similar procedure yields the expansion of the Dirac operator 



D{p,q) = 
Vip,q) = { 

+ { 
+ 



Do{p)6pg + V{p,q) 
1 



(39) 



(40) 



uj{p) + uj{q) 
1 



-Xl{p,q) 



}[MP,<l)~^Xl{p,q) 



Mp) 

a;(p) + a;(g)J L""^^''^^ ^(p) 



uj{q) 
Xoiq) 



uj{q) 



1 



;}{; 



;}{; 



■} 



X 



t^uj{p) + uj{q) ^ ^uj{p) + uj{t) J '^ooit) + uj{q) 
-Xo(p)X|(p,t)Xi(t,g) 

-Xi(p, t)Xo^(t)Xi(t, g) - Xi(p, t)Xlit, g)Xo(g) 
u;(p) + u;(t) + iu{q) 



u){p)uj{t)uj{q) 



-Xoip)Xlip,t)Xoit)Xlit,q)Xoiq) 1 + ■ ■ ■ . (41) 



This expansion may be obtained through the integral representation of the inverse 
square root of X^X: 



1 dt 1 



V^XtX J-oo 7f t^+X^X' 

Up to the second order, we have 

~ I ~To (-'^0-^2 + xIXq + XlXi 

J-oo n 4- X\Xn ^ 



(42) 



1 — (xlx-i + X}X()] 1 — (xlXi + XlXn] 1 — . 

^oont^ + XlXo^ ° U^ + XlXo^ ° h^ + X},Xo 

(43) 



Noting that (t^ + XqXq) ^ is diagonal in spinor space and it commutes with Xq, we 
obtain 



aD = aD, 







~l~ / 1 (t^X^ — XnX} Xn , 

J-oo TT t2 + X^o^O ^ " ^ t2 + X'oXo 

~\~ / 1 (t^Xo — XnXlXr)] 1 

J-oo nt^ + xlXo ^ ' ^ t^ + XlXo 

— ^ ~ Tprzr (^i) 1. TTTTT (^0^1 + ^1^0 



vr t2 + xlX^ 't^ + XlXo ^ ' + XlXo 

J-oo TT + xlXo ^ XlXo ' + Xo^Xo 

+ r ^ (Xo) ^ (xJXo) ^ fx|Xo) ^ + • ■ -(44) 

Going to the momentum space, the integration over the parameter t can be performed 
explicitly and we obtain Eq. (^Ol). 

The tree level propagator is given by 

Kp) = - E (1 - cos ap,) - -Mo. (46) 

CL ,, Qj 
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Since u^p) is positive definite, the pole of the propagator occurs when u{p) + b{p) = 0, 
which is fulfilled only if sin^ ap^ = and b{p) < 0. Therefore for — 2r < Mq < the 
propagator exhibits a massless pole only when p^ = 0. 

Now using the expansion Eq. ( PD we compute —atie'-f^D. (Here the local transfor- 
mation is treated so that e is space-dependent.) The non- vanishing contribution comes 
only from the term containing Xo{p)Xi{p,t)Xo(t)Xi(t, q)XQ{q) . We obtain 

-atie^^D = a' f e^e'^P-'>^^"'ghi{A^{p)A,{-q)}f^,{p,q), (47) 

Jpq 

ffiu{p,q) = h^{p + k,k,k + q) X 

J k 

ii^^Xo{k + p)d^Xl{k + p/2)X^{k)d,Xl{k + q/2)Xo{k + g), 

(48) 

nin k a) = { ^ ]{ ^ ]f ^ y uj{p)+uj{k)+uj{q) ^ 
^uj{p) + uj{q)^^u{p)+u{k)^^uj{k) + u{q)^^ uj{p)u{k)u{q) ^ 

(49) 

Then using the propagator S{p, ps) defined in five dimensional momentum space (p^, ps) G 
X R, 

S{p,p,)-' = Xoip)+t'y,p, (50) 

f^u{p,Q) is rewritten as 

f,,{p, q) = 2i[ ii{S{k + p, k^)d^S~\k + p/2, k,) 

J (ki,k'^)&T*xR 

S{k, k,)d,S-\k + g/2, A;5)5(A; + g, A;5)955-^(A;, ^5)}. (51) 

Now it is easy to verify the structure f^u{p,<l) = PpQTdpdrf^u{p,<l)\p,q=o + 0{a), and 
dpdrff^u{p,q)\p,q=o oc ep^^,,, leading to the anomaly g'^c{Mo,r)ep^rvFp^Fru when inserted 
in eq. ( ^Tj) with the coefficient 

c(Mo,r) = 2i^^epj^2p3^^^3 ^^^^trll/p^Lp^Lpgl/^^L^^j, (52) 

This expression shows that c(Mo,r) is invariant by a momentum-dependent, contin- 
uous change of the scale of the propagator S: S{p) Q{p)S{p) as long as fl{p) 7^ in 
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T"^ X R. It allows to replace S and ^ in eq. (0) with V and V ^ where 

V{p,p,) = N{p,p,)S{p,p5), N{p,p,) = ^pI + u{pY. (53) 

^ is a mapping from T"^ x i? to 5*^ and monotonously increases from —1 to 1 as 
increases from — oo to oo. This allows one to interpret that c(Mo,r) is the winding 
number of the mapping, V' , from to S*^ C derived from V by fixing p^ and it is an 
integer. For Mq > 0, riQ^cc b) is positive definite on x i? and the image of by V' 
does not cover C S^, leading to the value c(Mo, r) = 0. Further studies along this line 
leads to the value c(Mo,r) = l/lGvr^ for — 2r < Mq < |jl8[. (Note that our convension 
for the gamma matrices is as 7172737475 = — 1- ) 

This value is also obtained by evaluating the difference Ac = c(e, r) — c(— e, r) in 
the limit e — |]T^. A continuous and momentum dependent deformation of Mq does 
not change c{Mo,r) as long as 7^ in x i?, stemming from the fact of c(Mo,r) 
being the winding number. Here the surfaces Mq = ±e are deformed to m±{p) defined 
on T^, where m±{p) = for I]^=iP^ ~ and m±(p) = ±e for I]^=iP^ ~ S"^- Then the 
difference Ac arises from the integration in the vicinity of the center of the Brillouin zone 
in eq. (|52|). For sufficiently small 6 (still larger than e ), S and S^^ may be replaced by 
their continuum expressions, leading to 

by taking the limit e — for a finite small 6. 

Thus we obtain the final result of the axial U{1) anomaly as 

- atr {e^,D} = -£-^Nf j d'x e{x)e,,ruF;^{x)F!^,{x). (55) 

The same analysis can be done in two dimensions where the anomaly comes from the 
second term in eq. (^l]). The result is 

- atr {£75/}} = -^Nf j (fxe{x)e^,F^,{x). (56) 

We have discussed the weak coupling expansion of the lattice QCD with a Dirac 
operator satisfying the Ginsparg- Wilson relation using the expression of Ref. [^. We 
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confirmed that the anomalous behavior of the fermion partition function under the axial 
f/(l) transformation of Ref. ||^ is expressed in the form of the topological charge density 
for slowly varying perturbative gauge fields, which supplements the earlier calculation in 
Ref. i 

Narayanan, Vranas and Singleton Jr. have studied numerically the relation between 
the index of the Dirac operator Tr [H/^/lPj and the topological charge of lattice gauge 



field pO| . An interesting possibility for studying this type of Dirac operator in numerical 



simulation was considered recently by Chiu |^ . 
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